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TECHNICAL PAPER 


MORE ON EXACT STATE RECONSTRUCTION IN 
DETERMINISTIC DIGITAL CONTROL SYSTEMS 


I. INTRODUCTION 


Books on modern digital control systems usually address the problem of controlling a continuous- 
time plant driven by a zero-order-hold with a sampled output as shown in Figure 1. (For example, see 
Reference 1 .) A common solution to this problem is to reconstruct the state of the system at the sampling 
instants using a state observer and then feed back the reconstructed states [2], However, the state 
observer has two undesirable characteristics. First of all, it is a dynamical system in itself and, hence, 
adds additional states and eigenvalues to the system, which can affect system stability. Second, as a 
consequence, the reconstructed state is normally an approximation to the true state and is usually not a 
good one early in the state reconstruction process unless the initial state of the system is well known. 
Recently, Polites developed a new approach to state reconstruction which has neither of these problems 
[3], Subsequently, he extended this work and developed what he called the Ideal State Reconstructor 
[4,5]. It was so named because: if the plant parameters are known exactly, its output will exactly equal, 
not just approximate, the true state of the plant and accomplish this without any knowledge of the plant’s 
initial state. Besides this, it adds no new states or eigenvalues to the system. Nor does it affect the plant 
equation for the system in any way; it affects the measurement equation only. It is characterized by the 
fact that discrete measurements are generated every T/N seconds and input into a multi-input/multi- 
output moving-average (MA) process [6]. The output of this process is sampled every T seconds and 
utilized in reconstructing the state of the system. In this paper, a special form of the Ideal State Recon- 
structor is presented which is simpler to implement than the most general form. Before presenting this 
simpler form, some pertinent results to date for continuous-time plants driven by a zero-order hold are 
reviewed in Section II. Then, the Ideal State Reconstructor, in its most general form, is summarized in 
Section III. Finally, the special form of it is presented in Section IV. An example of this special form is 
given in Section V. The conclusions and some recommendations for future study are presented in Section 
VI. 



Figure 1 . Continuous-time plant driven by a zero-order-hold 
with standard measurements. 








II. PRELIMINARY 


For the plant in Figure 1 , x(t) e R n is the state vector, u(kT) e R r is the control input vector, ys(kT) 
€ R m is the standard output or measurement vector, F e R nxn is the system matrix, G e R nxr is the control 
matrix, and C s e R mxn is the standard output matrix. It is well known that this system can be modeled at 
the sampling instants by the discrete state equations [1] 


X(k + 1 ) = Ax(k) + Bu(k) 

(1) 

y s (k) = C s x(k) , 

(2) 


where 


4>(t) =c€~' [(si-F)' 1 ] , 

(3) 

A = 4>(T) , 

(4) 

B = / 4>(\) d k G . 

0 

(5) 


4>(t) e R nxn is the state transition matrix. A e R nxn is the system matrix and B e R nxr is the control matrix for 
the discrete state equations (1) and (2). 

A and B can be determined analytically using equations (3) to (5). An alternative approach, which 
is also quite suitable for numerical computation, is as follows [7]: 4>(t) and /' <J)(\)d\ can be expressed in 
the form of matrix exponential series as 0 


4>(t) 


3C 


V F‘ V 



(6) 


2 



and 


/ 4>(\)d\ = S 

0 i = 0 


f r' 
(i+l)! 


(7) 


respectively. From equations (6) and (7), 


<Kt) = I + F 



4>(X)dX 


(8) 


T 

where I e R nxn is the identity matrix. Hence, f 4>(\)d\ can be determined using equation (7) with t = T 
and this result substituted into equation (8) to get 4>(T). With these results, A and B can be found using 
equations (4) and (5). 


Now consider the plant in Figure 2, which is a generalization of the one in Figure 1 . In addition to 
the standard output y s (kT), the plant in Figure 2 has the output y F '(kT) generated as follows. First, the 
continuous-time output z(t) e R p is sampled every T/N seconds. Every N samples are multiplied by the 

weighting matrices Hj e R qxp , j = 0,1, ,N- 1 , and then summed to generate the output y F (kT) e R q , 

every T seconds. Functionally, this is equivalent to passing the discrete measurements generated every 
T/N seconds through a multi-input/multi-output MA process with coefficient matrices Hj, j = 0,1,...,N- 
1 , and then sampling the output of this process every T seconds to generate y F (kT). Then y F (kT) has E_ 
u[(k-l)T] subtracted from it, where E_ e R qxr , to produce the modified MA-prefiltered measurement vec- 
tor y F '(kT) e R q . This is catenated with ys(kT) to form the total measurement vector y-r(kT) e R m+q . In 
Figure 2, C F e R pxn since z(t) e R p and x(t) e R n . 

Previously, Polites [8] showed that when 


E_ = H(3 

where H e R qx(Np) is given by 


(9) 


H = [H 0 j H, 


Hn-,] 


( 10 ) 
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Figure 2. Continuous-time plant driven by a zero-order-hold with standard 
and modified MA-prefiltered measurements. 













and (3 e R (Np)xr is 


/ <KMd\ 


r -(T/N) 


P = 


/ <j)(X)d\ 


' -(N-1KT/N) 

C F I / <KMd\ 


( 11 ) 


the discrete state equations for the plant in Figure 2 become 


x(k+l) = Ax(k) + Bu(k) 


( 12 ) 


r Ys(k) 
y T (k) = — 

L y F '(k) 


;k) i r c s i 
'(ku Ld_J 


= — x(k) = Ct x(k) 


(13) 


where D_ e R qxn is given by 


D_ = Ha 


(14) 


and a e R (N P )xn is 


a = 


r c F 4>(Q) 

c f 4»(- t/n) 

L C F 4>[-(N-l)T/N] 


(15) 


5 



From equation (13), 


C T 



( 16 ) 


where C T e R (m -‘i )xn . 

E_ and D_ can be evaluated analytically using equations (3), (9) to (11), (14), and (15). An 
alternative approach, which can be either analytical or numerical, is as follows. Let t = -j(T/N), where j 
= 0,1,...,N-1, and use equation (7) to determine J' j(T/N) 4>(X)dX, j = 0,1,...,N-1. Substitute these 
results into equation (8) to get 4>[-j(T/N)] , j = 0,1,..., N-l. At this point, E_ and D_ can be found using 
equations (9) to (11), (14), and (15). 


III. THE IDEAL STATE RECONSTRUCTOR 


A general block diagram of the plant and the Ideal State Reconstructor, in its most general form, 
is shown in Figure 3. Notice the similarity between Figures 2 and 3. By virtue of this, if E_ is given by 
equation (9), then equations (12) to (16) define the discrete state equations for the system in Figure 3 up 
to the output y T (k). Proceeding further, y T '(k) is related to y T (k) by the expression 


y T '(k) = (C t t C t )-' C T T y T (k) 


(17) 


However, for equation (17) to be meaningful, (C T T C T ) _I must exist, and this occurs only when (C T T C T ) 
is nonsingular. Recall that C T e R< m+( i> xn . If (m + q) 5= n and C T has maximal rank (i.e. , rank n), then (C T T 
C T ) is positive definite and therefore nonsingular [9]. Hence, equation (17) requires that (m + q) 2= n and 
rank (C T ) = n for it and the Ideal State Reconstructor to be meaningful. Assuming this is the case, it 
follows from equations (12), (13), and (17) that the discrete state equations for the system in Figure 3 are 


x(k + 1 ) = Ax(k) + Bu(k) 

(18) 

y T '(k) = x(k) • 

(19) 


Hence, the output of the Ideal State Reconstructor, y T '(kT), exactly equals the true state of the plant, 
x(kT). Consequently, if one is given the plant in Figure 1 and modifies it to conform to Figure 3, he can 
exactly reconstruct the state of the plant without adding any new states, eigenvalues, or dynamics to it, 
since the plant equation (18) for the system in Figure 3 is identical to the plant equation (1) for the plant in 
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Figure 3. Block diagram of the 


IDEAL STATE RECONSTRUCTOR 


1 





and the Ideal State Reconstructor in its most general form. 











Figure 1 . In Figure 3, exact state reconstruction is achieved when E_ is given by equation (9) and C T is 
given by equation (16) where D_ is given by equation (14). In addition, C T must satisfy the requirements 
just imposed on it. 

One of these is that the dimensions of C T , namely (m + q)xn, satisfy the relationship (m + q) > n. 
However, this can be rewritten as q 3= (n-m). Hence, the number of rows, q, in the weighting matrices 
Hj, j = 0,1 ,...,N-1 , must equal or exceed the number of states, n, in the plant state vector, x.(t) or x(kT), 
minus the number of elements, m, in the standard measurement vector, y s (kT). Since q can be chosen 
arbitrarily, this requirement can be easily satisfied. The other requirement on C T is that it have rank n. 
One approach to satisfying this requirement is as follows. Recall that C T is defined by equation (16) 
where D_ is defined by equation (14). Assuming C s is given, then one can choose D_ so that C T has rank 
n and then find H to give the desired D_. One solution to the problem of finding H to give the desired D_ 
is to let 


H = D_ (a T a )' 1 cx T . 


( 20 ) 


This follows from equation (14). However, like before, this requires that (a T a) be nonsingular. Recall 
that a e R (Np)xn . if (Np) 3= n, or equivalently N 2= n/p, and a has maximal rank (i.e. , rank n), then (a T a) 
is nonsingular. The first requirement can be easily satisfied because the number of weighting matrices 
N, where the weighting matrices are Hj, j = 0,1,..., N-l, can be arbitrarily chosen so that N ^ n/p. 
Recall that n is the number of states in the plant state vector, x.(t), and p the number of elements in the 
output vector ^(t). 

In summary, the procedure to achieve exact state reconstruction with the Ideal State Reconstruc- 
tor is as follows. Given the plant in Figure 1, modify it to conform to Figure 3. Choose the number of 
rows, q, in the weighting matrices Hj, j = 0,1 ,...,N-1, so that q ^ (n-m) where n is the number of states 
in the plant state vector, x(t) or x(kT), and m is the number of elements in the standard measurement 
vector, y s (kT). Choose the number of weighting matrices, N, so that N 3= n/p where p is the number of 
elements in the output vector z(t). Assuming the (Np)xn matrix a, defined by equation (15), has maximal 
rank (i.e., rank n), let H be given by equation (20) where D_ is chosen so that the (m + q)xn matrix C T , 
defined by equation (16), has maximal rank (i.e., rank n). The weighting matrices Hj, j = 0,1,...,N-1, 
are found by partitioning H as in equation (10). Finally, let E_ in Figure 3 be given by equation (9) where 
(3 is defined by equation (11). The discrete state equations for the system in Figure 3 are now given by 
equations (18) and (19). Hence, the output of the Ideal State Reconstructor, y T '(kT), exactly equals the 
true state of the system, x(kT). 

In the event no standard measurements are used in the state reconstruction process, then_y s (kT) is 
a null vector, C s is a null matrix, and the Ideal State Reconstructor in Figure 3 degenerates to the one 
presented in Reference 4. In this case, the methods in Reference 4 for choosing the parameters in the 
Ideal State Reconstructor to achieve exact state reconstruction, as well as the one described here, are 
applicable. 


8 



IV. A SPECIAL FORM OF THE IDEAL STATE RECONSTRUCTOR 


As indicated in Section III, one requirement of the Ideal State Reconstructor is: the number of 
rows, q, in the weighting matrices Hj, j = 0,1 must be chosen so that q 2= (n-m) where n is the 

number of states in the plant state vector, x(t) or x(kT), and m is the number of elements in the standard 
measurement vector y s (kT). A special form of the Ideal State Reconstructor, which is simpler to 
implement than the most general form, can be obtained by letting q = n-m. In general, C T e R""' 1 ' 1 '" . 
Hence, for this special case, C T e R nxn , which is a square matrix. Now, if D_ can be chosen so that C T is 
nonsingular, then 


y T '(k) = Ct 1 y-r(k) 


( 21 ) 


which follows from equation (17), and the Ideal State Reconstructor assumes the special form shown in 
Figure 4. From equations (12), (13), and (21), the discrete state equations for the system in Figure 4 
become equations (18) and (19), like before. Again, the output of the Ideal State Reconstructor, y T '(kT), 
exactly equals the true state of the plant, x(kT), but with a minimum number of rows in the weighting 
matrices, in this case. 

Observe, from equations (12) and (13), that when Ct is square and nonsingular, y T (k) becomes a 
legitimate state vector for the system. Solving for_x(k) in equation (13) and substituting this result into 
equation (12) yields the plant equation for the system in terms of the state vector y T (k). The result is 


y T (k+l) = A* yr(k) + B* u(k) 


where A* = C T A C T _I and B* = C T B. Consequently, if it is sufficient to construct the state vector 
y(kT) in place of the state vector x.(kT), then the special form of the Ideal State Reconstructor in Figure 4 
can be simplified, because there is no reason to calculate y T '(kT). 


V. AN EXAMPLE 


Consider the double integrator plant and the special form of the Ideal State Reconstructor shown 
in Figure 4. Manipulating this plant into the format of Figure 3 yields 


F 



1 

0 


( 22 ) 
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PLANT 


IDEAL STATE RECONSTRUCTOR 



I -1 


Figure 4. Block diagram of the plant and the special form of 
the Ideal State Reconstructor. 
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Figure 5. Double integrator plant and the special form of the 
Ideal State Reconstructor. 










G 


0 


(23) 


and 

C s = C F = [1 0] . (24) 

Since F e R nxn , G e R nxr , C s e R mxn , and C F e R pxn , it follows from equations (22) to (24) that n = 2 and r 
= m = p = 1. Using equations (22) and (23) and the formulas presented in Section II, 



(25) 


/ cJ>(\)dA. 

o 



(26) 



T 2 /2 

B = 

T 

For the special form of the Ideal State Reconstructor, q = n-m = 1 . The requirement N 5= n/p can be 
satisfied by letting N = 4. Now a and (3 can be evaluated using equations (11), (15), (23) to (25), and 
(26), and are found to be 


12 



0 


1 


a 


1 

1 

1 


T 

4 

2T 

T 

3T 

4 


and 


0 = 


0 

T 2 

32 

4T 2 

12 

9T 2 

32 


1 


(27) 


(28) 


respectively. In equation (27), eliminating any two rows forms a 2x2 matrix with nonzero determinant, 
assuming of course T > 0. Hence, rank (a) = 2 = n and so (a r a) is nonsingular. Consequently, 
(a T a)'' a T exists and is found to be 


(a T a)' 1 a T 



('.) (lo) 

(i) ( 4 ) 



(29) 


using equation (27). Since (a T a)" 1 a T exists, H can be given by equation (20) where D_ needs to be 
chosen so C T is square and nonsingular. Since D_ e R qxn , then D_ becomes a 1x2 matrix. Now, if D_ is 
chosen to be 


D_ = [0 1] 


(30) 


it follows from equations (16), (24), and (30) that 


C T = C T -‘ 


1 0 
0 1 
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Obviously, C T is square and nonsingular and there is no need to calculate y T '(kT) in Figure 5. From 
equations (10), (20), (29), and (30), 



/ 6 \ 

/ 2 \ 

, ( 2 \ 

r 

H = [H 0 j H, j H 2 1 H 3 ] = 

\5t) 

( 5T ) 

! \ 5T / 

' V 5T / j 


which reveals the weighting matrices Hj, j = 0,1, 2, 3. From equations (9), (28), and (31), 


E_ = -3T/8 


The special form of the Ideal State Reconstructor is now completely defined for this example. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


This paper has presented a special form of the Ideal State Reconstructor for deterministic digital 
control systems which is simpler to implement than the most general form. The Ideal State Reconstructor 
is so named because: if the plant parameters are known exactly, its output will exactly equal, not just 
approximate, the true state of the plant and accomplish this without any knowledge of the plant’s initial 
state. Besides this, it adds no new states or eigenvalues to the system. Nor does it affect the plant equa- 
tion for the system in any way; it affects the measurement equation only. It is characterized by the fact 
that discrete measurements are generated every T/N seconds and input into a multi-input/multi-output 
MA process. The output of this process is sampled every T seconds and utilized in reconstructing the 
state of the system. 

The Ideal State Reconstructor is ideally suited for systems where measurements are available at a 
faster rate than the control law equations need to be solved. A good implementation of it would be to 
have a microprocessor dedicated to solving the MA-prefilter calculations recursively as the 
measurements become available every T/N seconds. Every T seconds or N calculations, the result could 
be transferred to a central processor where the remaining calculations in the reconstructor are made and 
the control law equations are solved. 

Since it is an open loop type state reconstructor, the Ideal State Reconstructor may be less robust 
than the state observer when parameter uncertainties and measurement and process noise are considered. 
This is an area for further study. However, intuitively it would seem: the more measurements used in 
reconstructing a given state, the more robust the reconstructor. If so, then making N as large as possible 
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would help. Perhaps, N could be as large as 100, or even 1000, with the suggested implementation, in 
many practical problems. If it turns out that robustness is still a problem, then the future of the Ideal State 
Reconstructor may be that of a prefilter to the state observer, or better yet the Kalman filter. In a system 
where measurements are available at a much faster rate than the Kalman filter algorithms can be solved, 
the Ideal State Reconstructor might provide a descent estimate of the state which could then be further 
improved by the Kalman filter. This also is an area for future study. 
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